ON THE COMPLEMENTARITY OF THE QUADRATURE OBSERVABLES 
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Abstract. In this paper we investigate the coupling properties of pairs of quadrature observ- 
ables, showing that, apart from the Weyl relation, they share the same coupling properties 
as the position-momentum pair. In particular, they are complementary. We determine the 
marginal observables of a covariant phase space observable with respect to an arbitrary rotated 
reference frame, and observe that these marginal observables are unsharp quadrature observ- 
ables. The related distributions constitute the Radon tranform of a phase space distribution 
\ of the covariant phase space observable. Since the quadrature distributions are the Radon 

transform of the Wigner function of a state, we also exhibit the relation between the quadra- 
\ ture observables and the tomography observable, and show how to construct the phase space 

£^ ■ observable from the quadrature observables. Finally, we give a method to measure together 

with a single measurement scheme any complementary pair of quadrature observables. 
PACS number: 03.65-w,03.65.Ta,0365.Wj 
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Dedicated to Peter Mittelstaedt in honour of his eightieth birthday. 



d • 1. Introduction 

3 . 

The notion of complementarity was introduced to the vocabulary of quantum physics by Niels 
Bohr in his famous Como lecture of 1927 as a key to the understanding of quantum phenomena 
in terms of classical concepts [5]. One of the most explicit uses of the "traffic rules" of Bohr 
was in his 1935 paper [6], where he argued that position and momentum of a particle are 



complementary quantities in the sense that all the experimental arrangements allowing their 
unambiguous operational definitions are mutually exclusive but they both are needed for a full 
description of the situation. 

In addition to the position-momentum pair, energy-time, path-interference, number (action)- 
phase, spin-phase, or spin components, are frequently occuring examples of pairs of complemen- 
tary observables. Moreover, complementary modes of description, like the use of complementary 
bases, or the past and the future state determinations of the system, are often discussed cases. 
^ ■ For an overview of various aspcets of this notion we refer to [TOl 18]. 

In this paper we investigate the properties of pairs of quadrature observables in close analogy 
to the position-momentum case, reviewed in section [2l We show in section [3] that - apart 
from the Weyl relation - they share all the coupling propeties of position and momentum, 
which reflect the strong incompatibility of these observables. In sections [4] and [5] we exhibit 
the connection of the quadrature observables to the covariant phase space observables and 
the tomography observable, respectively. In the concluding section [6] we demonstrate, that 
though any pair of quadrature observables is complementary in the sense that none of their 
measurements can be combined into a joint measurement of theirs, there are single measurement 
schemes which allow one to determine the measurement outcome distributions of the given 
quadrature pair for a large class of states of the system. 

2. The pair (Q, P) 

Let 7i = L 2 (R) be the usual L 2 -function space on R spanned by Hermite functions h n , 
n G N = {0,1,...}. Consider the selfadjoint position operator Q on the Hilbert space TC of 
a quantum object in one dimension, and let Q be its spectral measure, so that Q(X) is the 
multiplication with the characteristic function xx of the (Borel) set X G £>(R). Let F be the 
unitary Fourier-Plancherel operator on L 2 (R) so that P = F~ l QF is the selfadjoint momentum 
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operator —id/dx, with the spectral measure P = F 1 QF. For any state p (positive trace-1 
operator) we let p Q and p p denote the densities of the probability measures X i— ► P®(X) = 
tr [pQ(X)] and Y i— > Pp(X) = tr [pP(Y)] with respect to the Lebesgue measure. 

The unitary groups {e~ igP | q G M} and {e ipQ | p G M} of Q and P fulfill the Weyl relation 

m e -iqP e ipQ — e -iqp e ipQ e -iqP 

and, modulo unitary equivalence, the pair (Q,P) is uniquely determined by this relation [27] . 
In addition, the pair (Q, P) has the following well-known coupling properties: 

a) QP — PQ = il (on the (dense) domain of the commutator); 

b) inf {Var(p Q ) • Var(p p ) | p a state } = \ > 0; 

c) com(Q,P) = {ip(E L 2 (R) \ Q(X)P(Y)V> = P(Y)Q(X)ip for all X, Y G B(R)} = {0}; 

d) Q(X) A P(Y) = Q(X) A P(R \ Y) = Q(R \ X) A P(Y) = for all bounded X, y G B(R); 

e) tr[Q(X)P(Y)] = ^A(X)A(Y) for all bounded X, Y G with A(X) being the 
Lebesgue measure of X. 

In addition to the Weyl relation, all the properties a) through e) reflect extreme incompatibility 
of position and momentum observables. The commutation relation a) as well as the preparation 
uncertainty relation b) belong to the basic arsenal of quantum mechanics and need no further 
comments here. The property c) expresses the fact that for no state p the map (X, Y) \— > 
tr [pQ(X) A P(Y)] extends to a probability measure on £>(M 2 ) [29j. On the other hand, relations 
d) and e) have been taken to describe the complementarity of these observables in the sense 
of lack of any joint measurements, see, for instance, [HI III. 8. .2], [8, IV. 2. 3], or, as a kind 
of generalization of the "complementary bases" of the finite dimensional case [U [20]. It is, 
perhaps, well-known, and will also be shown below that none of these five properties a) — e) is 
sufficient to determine the pair (Q, P) to be the Weyl pair, for an explicit proof of the case e), 
see, e.g. [T4] . 

It is also well-known that the pair (Q, P) is informationally incomplete: the measurement 
outcome statistics p Q ,p p of these observables do not suffice, in general, to determine the state 
p of the systemQ In the words of C.F. von Weizacker [28], this is a reflection of the surplus of 
information coded in the quantum notion of state, when compared with the classical one. 

3. The pair (Q, Q e ) 

Position Q and momentum P can be obtained in a smooth way from each other. Indeed, let 
Ug = e l6H , 6 G R, be the unitary operator defined by the oscillator operator H = \{Q 2 + P 2 ), 
and define Qg = UgQUg, so that Qe is the quadrature operator, with the spectral measure 
Qe = UgQUg. Clearly Q = Q and = P (since F = in fact, Qg = Qcos9 + Psin9. 

To study the coupling properties of any two quadratures (QanQp), it is sufficient to consider 
the pair (Q,Qe), G [0,27r), since for any pair (Qa,Qp) one finds a unitary operator U := U a 
such that Q a = UQU* and Qp = UQ/3- a U*, that is, the pair (Q a ,Qp) is unitarily equivalent 
to the pair (Q, Qg) where 6 = (3 — a. 

Using the operator relation 

giy(Q cos 8+P sin 9) ^iyQ cos O^iyP sinO ^iy 2 (cos 6 sin 9)/2 

together with the Weyl relation one checks that the unitary operators e lx ® and e %yCSe fulfill the 
Weyl relation exactly when 9 = n/2, that is, Qg = P. 
Clearly, for any pair of quadratures (Q,Q$), 

QQe - QeQ = isin9 1 



One of the first examples demonstrating this fact is reported in [26] and is due to V. Bargmann. 
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on the dense domain of the commutator. Denote b := sin 9 and assume that b ^ 0. Since 
(tQ, Qe) is a Weyl pair, and the spectral projections of are of the form Q(bX), X G S(R), 
one notes that all the coupling properties b) — e) hold for the pair (Q, Qe), as well. In particular, 
for any 9 G" {0, 7r}, 

sin 2 9 

(2) inf {Var(p Q ) ■ Var(p Qfl ) | p a state } = — ^— > 0, 

the quadratures Q and Qe are totally noncommutativdl, they do have no joint probability of 
the form (X, Y) \— > tr [pQ(X) A Qg(Y~)], they are complementary in the sense of d) and they 
satisfy the trace formula e). 

As already pointed out, the pair (Q, P) is informationally incomplete. Clearly, the same is 
true for any pair (Q,Qe). These pairs can, however, be completed adding further quadratures. 
Indeed, any set {Qe \ 9 G S}, where S C [0, 7r) is a dense set, is informationally complete, and 
thus allows state determination on the basis of the statistics p Qs , 9 G S, see, e.g. [T3l [T8] . The 
other method to complete the pair (Q, P) is to replace it by a coexistent pair (p * Q, v * P) of 
unsharp position and momentum observables such that their joint observable is informationally 
complete [3] 0]. We recall that, for instance, p * Q is the normalized positive operator measure 
(POM) defined by the convolution of the probability measures p and , 

tr[p(p*Q)(X)] = (p*pQ)(X)= / ^{X-q)dp Q p {q) = [ p(X - q)p Q (q)dq 

Jr Jr 

where p is a state. As will be seen below, these two approaches are closely related with each 
other. 

4. The pairs (Q, Q 9 ) and the phase space pom 

Due to the noncommutativity of the pair (Q,Qe), there is no normalized positive operator 
measure (POM) E : £>(IR 2 ) — > C{7i) which would have both Q and Qe as the marginal observables 
[24, Thm IV. 1.3.1.]. On the other hand, the unsharp pair (p * Q, v * P) has joint observable 
exactly when the probability measures p and v have Fourier related densities, in which case 
(p * Q, v * P) are the (Cartesian) marginal observables of a covariant phase space observable 
Gk generated by a positive trace- 1 operator K on 7i [TJ]. We recall that Gk is defined by the 
operator density (q,p) i— > W (q, p)KW (q, p)* , that is, 

G K {Z) = J- / W(q,p)KW(q, P ydqdp, Z G B(R 2 ), 
J z 

where W(q,p) = e tqp ^ 2 e~ tqP e ipQ is the Weyl operator. The Cartesian marginal observables 
X i — > G K (X x R) and Y i— > G K (R x Y) are, indeed, of the form p K * Q and u K * P, with the 
convolving probability measures p K {X) = p® KIl *(X) and u K (Y) = PyikwO^)^ where II = II* = 
U± n is the parity operator (11-0) (x) = x). For any (bounded) operator A we will use the 
following short notation: 

A e := U e AU;. 

Especially, HKXl* = K n . Consistently with the notation p Q , the densities of p K and v K may 
be written as = {UKU*) Q = K n * Qn = K Q ~ and = K p * . 

One may also determine the marginal observables of Gk with respect to a rotated orthonormal 
frame {ei(9), e 2 (6 1 )} of R 2 where 

e^^) := (cos#, sin^), e 2 (#) := (— sin 9, cos 9) 

so that qe^i(9) + pee 2 (9) = (q,p) with 

qe = q cos 9 + p sin 9, pe = —qsm9+p cos 9. 



2 We recall from [16] that condition J5]) alone implies that com(Q,Qg) = {0}. 
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Since UgW(q,p)Ug = W(qcos9 + psm9,—qsm9 + pcos9) = W(qg,pg), the relevant marginal 
observables are simply 

X _L f W(q, P )KW(q, P yd Pe = (fi K -° * Q e )(X), 

Z7V JxxR 

y^l^l W(q,p)KW(q,p)*dq e = (v K -° * P e )(Y), 

Z7T JrxY 

that is, they are unsharp quadrature observables, the convolving measures being determined 
by the rotated generating operator K_g. For any state p, the density of \i K ~ B * Qg, say, is the 
convolution K®_ e * p Qe of the densities K%_ $ = (UK_ e U*) Q and p Qe = p Q e . 

The "^-marginal observables" p Ke * Qg of constitute the Radon transform Gk- Indeed, 
for any state p, let g K (q,p) '■= tr \pW(q,p)KW(q,p)*] denote the density of Gk in the state p. 
The Radon transform of g p K is defined as 

(Rg p K )(9, q e ) : = / g p K (q e cos9 - p e sin 9, q e sin 9 + p e cos 9) dp e 
Jr v v ' 

= (g,p) 

This shows that for any 9 G [0, 27r), the function qg \— > (Rg p K )(9,qg) is (2tt times) the density 
K®_ e * p Qs of the probability measure p K -° * p® B of the unsharp rotated quadrature observable 
p K - * Qg in the state p. 

Recall that, for fixed 9 and qg, the image 1(9, qg) of 

K 3 pg i— > (qe cos 9 — pg sin 0, sin 9 + pg cos 0) G M 2 

is a line on the plane M. 2 . It goes through a point (qg cos 9 , qg sin 9) = qgei(9) and its direction 
unit vector is e 2 (9). Hence, in the above definition of R, the integral is a line integral over 
1(9, qg) = q$ei(9) + IRe 2 (#) which is perpendicular to the vector qge\(9). 
Fixing 9, let X G B(R), and define 

Z(9,X) := { qeei (9)+pge 2 (9) G R 2 | qg G X, p e G M} = |J £(9,q e ) G B(R 2 ). 

qe ex 

Then 

ti[pG K (Z(9,X))} = ^ / g P K (q,p)dqdp=±- [ (Rg P K )(9 , qg)dq e 

to JZ(6,X) to J x 

= [ (K^g*p Q °)(qe)dqg. 
Jx 

For example, in the case of the number state K = \h n )(h n \, one gets 

G ]hn){hnl (Z(9,X)) = — i-= / / [ff n (a; - g ,)]V<*-»> a Q,(<fc)d© 

(where H n is the nth Hermite polynomial) and especially 

G lho){hol (Z(9,X)) = ^= [ [ e-^ 2 Qg(dx)dq e . 
V n Jx Jr 

Let W p be the Wigner function of the state p, that is, \N p (q,p) = Hx [pW(q,p)UW(q,p)*]. As 
well-known, the Radon transform of the (integrable) Wigner function of a state is the rotated 
quadrature distribution of this state, that is, (R\N P )(9, x) = p Qe (x), see, for instance, [23j. 
Therefore, the density x i— > (Rg p K )(9,x) is a smearing of the density x \— > (R\N p )(9,x) with the 

density K®_ 9 . This observation bring us to the tomography POM E ht . 
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5. The pairs (Q, Q e ) and the tomography pom E ht 

For any state p, the random sampling of the distributions p Qe , 6 G [0, 27r), determines a 
probability bimeasure (Q,X) h- > f Q f x p Qf, (x)^|^, which, when taken all together, determine 
the tomography observable, 

E ht (0 xl):=i/ Q fl (X)d0, 

studied extensively, for instance, in [2]. In particular, E ht is informationally complete. Its 
marginal observables are i— *■ ^ J Q dd I and 

r^e(X)de = J2 I [h n (x)] 2 dx\h n )(h n \ 

^ ^° n=0 ^ X 

since Q (X) = En,m=o ^" )d J x K{x)h m {x)dx \h n )(h m \. 

The statistics of a phase space observable can be obtained from the statistics of the 
tomography observable Eht in terms of a generalized Markov kernel, at least whenever the 
generating operator K is smooth, that is, the integral kernel of K belongs to the Schwartz 
space of R 2 . Indeed, in that case one may write for any compact Z C R 2 

-27T 



G K (Z)= f f A- f M« p (e,x)dqdp 
Jo Jr L Z7r Jz 



dE ht (9,x) 



where M^ p (9,x) is a smooth function with respect to all variables (25} sect. 3.2]. (We call 
this function a generalized Markov kernel since it is not necessarily positive.) For instance, if 
K = \h n )(h n \, the function M^ p (6,x) takes the form (25] eq. (3.7), Thm 1] 

M^>^(9,x) = M h n){hnl (0,x-q e ) 



,A fn\ 2 1 -" d 2u+1 [ _ {x _ qe) 2 f x - qe y2 
ul 9x2u+1 l Jo 6 



u=0 

oo 



2^\n) 2«(2k)\ H2k{x 



k=n 



For any smooth trace-class operators K and p we have [25], eq. (3.3)] 



so that 



(K^_ e *p Q o)(q 6 )= [ K Q °(x-q 6 )p Q °(x)dx = -±-- [ [ [ M« p (6' ,x)p Q °' (x)d6'dxd P g. 
Jr k^) Jr Jr Jo 



n\ ? 



The phase space observable G^, with K commuting with N, that is, K = Y^=o w nP[h 
< w n < 1, Y^=o w n = 1, is of special interest since then the angle marginal observable (with 
respect to the polar coordinates) of Gk is phase shift covariant, that is, a phase observable [21, 
thm 4.1]. Clearly, in that case Kf = K Q for all 6 G R. If K is smooth and commutes with N, 
then M ^ (#,x) = Mg o (0,x) for all x G R. 

6. Measuring the pairs (Q, Qg) 

As already pointed out, there is no POM having Q and Qg as its marginal observables. These 
observables do not have any joint measurements. Apart from that there are single measurement 
schemes which allow one to determine both the Q and the Qg -distributions p Q and p Qe for a 
large class of states p. 

To illustrate this possibility, consider a sequential combination of the standard von Neumann 
measurements of first Q and then Qe, as described, e.g. in [TT], III. 2. 6]. Such a sequential 
measurement defines a unique phase space observable G [IS], [7]. Its first marginal observable 
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is an unsharp position p * Q defined by the first measurement, whereas its second marginal 
observable is an unsharp quadrature v * Qg defined by the second measurement under the 
influence of the first measurement. The structure of the convolving measures p and v depend 
on the details of the applied measurement schemes, in particular, of the initial states of the 
probe systems. We do not need these details here. 

For any state p one may determine the moments of the marginal distributions p * p® and 
v * p® 6 , of the actual measurement statistics p^, and they are of the generic form 

(A* **?)[*]= / Jdfa*p*)(x) = ^(%[k-n]p*[n], 

where, for instance, p,[k) denotes the k th moment of p. Choosing the initial states of the two 
probe systems such that all the moments of p and v are finite (for instance, choosing the two 
states to be Gaussians) and assuming that p = \ip)(ip\, with ip in the linear hull of the Hermite 
functions, then the actually measured moment sequencies {{fJ>*p < j})[k])ken an d {iy *Pp 9 )[^])fceN 
can be solved for the sequencies (p®[k])km an d (p® e [k])k€N, which, due to exponential bound- 
edness of the involved probability measures, uniquely determine the distributions p® and p® 6 , 
respectively; for technical details, see 0,[l9]. Note that the distributions p Q and p Qe , with the 
above choice of p, suffice to determine the whole observables Q and Qg, respectively. Note also 
that the phase space observable G is not of the form unless 9 — | . 

To close this section, we recall that the eight-port homodyne detector with a strong local 
oscillator is an actual quantum optical implementation of a single measurement scheme which 
allows one to determine the distributions p Q and p Qe of any pair of quadratures (Q, Qg) for a 
large class of states [T7] . 

7. Concluding remarks 

We have studied the pairs of quadrature observables (Q, Qg), showing, in particular, that they 
share all the familiar coupling properties of the position-momentum pair (Q, P), except their 
defining property of being a Weyl pair. We have also determined the 6*-marginal observables of 
a covariant phase space observable Gk, and they turned out to be unharp quadrature observ- 
ables, the convolving probability measure being determined by the rotated generating operator 
K_g. These marginal observables p K ~ 6 *Qg constitute the Radon transform of the phase space 
observable Gk- Since the Radon transform RW p of the Wigner function W p of a state p gives 
the quadrature distributions p Qe in that state, we also exhibited the construction of the phase 
space observable Gk in terms of the tomography observable E^t defined as a random sampling 
of the quadrature distributions p Qe , 9 G [0,2n). We also showed that in spite of the fact that 
the quadrature observables Q and Qg are complementary observables in the sense that they 
have no joint measurements, it is, anyway, possible to measure the two observables together 
with a single measurement scheme, using, for instance, the statistical method of moments. 
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